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ABSTRACT 

Characterizing the nature and spatial distribution of the lensing objects that 
produce the previously measured microlensing optical depth toward the Large 
Magellanic Cloud (LMC) remains an open problem. We present an appraisal of 
the ability of the SuperMACHO Project, a next-generation microlensing survey 
directed toward the LMC, to discriminate between various proposed lensing pop- 
ulations. We consider two scenarios: lensing by a uniform foreground screen of 
objects and self-lensing by LMC stars. The optical depth for "screen-lensing" 
is essentially constant across the face of the LMC, whereas the optical depth 
for self-lensing shows a strong spatial dependence. We have carried out exten- 
sive simulations, based upon data obtained during the first year of the project, 
to assess the SuperMACHO survey's ability to discriminate between these two 
scenarios. In our simulations we predict the expected number of observed mi- 
crolensing events for various LMC models for each of our fields by adding artificial 
stars to the images and estimating the spatial and temporal efficiency of detect- 
ing microlensing events using Monte-Carlo methods. We find that the event rate 
itself shows significant sensitivity to the choice of the LMC luminosity function, 
limiting the conclusions which can be drawn from the absolute rate. If instead 
we determine the differential event rate across the LMC, we will decrease the 
impact of these systematic biases and render our conclusions more robust. With 
this approach the SuperMACHO Project should be able to distinguish between 
the two categories of lens populations. This will provide important constraints 
on the nature of the lensing objects and their contributions to the Galactic dark 
matter halo. 

Subject headings: dark matter — microlensing — galaxies: structure — galaxies: halos — 
Magellanic Clouds — galaxies: individual (MWG) 
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1. Introduction 

An elegant way to further our understanding of dark matter halos and to search for 
astrophysical dark matter candidates is to utiHze the defining feature of the dark matter: the 
effect of its gravitational field. Paczynski (1986) suggested searching for dark matter (DM) in 
the form of MAssive Compact Halo Objects (MACHO's) using gravitational microlensing. 
Several groups followed this suggestion and established microlensing searches toward the 
Large Magellanic Cloud (LMC) and other nearby galaxies. The MACHO group reported 
13-17 microlensing events toward the LMC (Alcock et al. 2000) with event timescales ranging 
between 34 and 230 days. They estimated the microlensing optical depth toward the LMC 
to be r = 1.21q 3 x 10~^. If we assume that MACHO's are responsible for this optical 
depth, then a typical halo model allows for a MACHO halo fraction of 20% (95% confidence 
interval of 8%-50%) with MACHO masses ranging between 0.15 and 0.9 M©. The OGLE 
collaboration also reported one LMC microlensing event (Udalski et al. 1997). Notably, none 
of the surveys toward the LMC have detected events with timescales Ihr <i < lOdays. This 
lack of short timescale events puts a strong upper limit on the abundance of low-mass DM 
objects: objects with masses lO^^M© < m < 10~^Mq make up less than 25% of the halo 
dark matter. Further, less than 10% of a standard spherical halo is made of MACHO's in 
the 3.5 X 10~ ''Mq < m < 4.5 x IO^^Mq mass range (Alcock et al. 1998). It is interesting to 
note that early results from surveys toward M31 are confirming a non- negligible MACHO 
content in M31's halo (Uglesich et al. 2004; Calchi Novati et al. 2005; Jetzer et al. 2004). 

Recent data and publications have reinvigorated discussion concerning microlensing 
rates toward the LMC. The EROS-2 project has found evidence for variability in an event 
classified as microlensing by the MACHO project (Tisserand & Milsztajn 2005) while also re- 
porting 4 new microlensing candidates toward the LMC. Additionally, using neural networks 
to analyze a subset of MACHO lightcurves including all events classified as microlensing by 
MACHO, Evans & Belokurov (2004) find an optical depth toward the LMC more consis- 
tent with that expected from known stellar populations. Bennett et al. (2005), however, 
show that previously unpublished data suppport the microlensing interpretation questioned 
by Evans and Belokurov. Accounting for new evidence that removes some MACHO events 
from the microlensing set, Bennett (2005) recalculated the microlensing optical depth toward 
the LMC using efficiencies determined for the entire data set. His revised optical depth is 
r = 1.0±0.3 X 10"^. He found the MACHO data to be consistent with a ~16% MACHO 
halo. EROS-2 calculate a preliminary optical depth of about 1.5 x 10^^ (Jetzer et al. 2004). 

Despite these varying constraints on the MACHO halo fraction, the microlensing event 
rate (reported by Alcock et al. 2000) toward the LMC significantly exceeds that expected 
from known visible components of our Galaxy. This raises the question of where the lenses 
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reside. Unfortunately the main observable in any given microlensing event, its duration, 
depends upon a combination of lens mass, position, and velocity relative to the source. Any 
conclusion about the spatial location of the lens population therefore depends upon the 
assumptions made about its mass and velocity. We note that in cases where the lightcurve 
exhibits a departure from the point-source/point-lens event shape (due to a binary lens, a 
binary source, non-inertial motion in the lensing system, etc.), this degeneracy can be lifted. 

Using the standard Galactic and LMC model, an optical depth of about 10^^ is expected 
toward the LMC from known Galaxy (e.g. thick disk, halo) and LMC components. This 
is significantly lower than the optical depth of 1.2 x 10^^ detected by the MACHO survey 
(Alcock et al. 2000). Given the difficulty associated with locating the lenses along the line of 
sight, previous microlensing surveys of the LMC have been unable to discriminate between 
a variety of possible sources for the excess LMC event rate. These include: (1) lensing by a 
population of MACHO's in the Galactic Halo, (2) lensing by a previously undetected thick 
disk component of our Galaxy, (3) disk-bar or bar-bar self-lensing of the LMC, or (4) lensing 
by an intervening dwarf galaxy or tidal tail. 

See § 2 for a more detailed discussion of these populations. Due to the limited number 
of events observed to date it is not yet clear which scenario or combination of scenarios 
explains the observed lensing signal. 

The SuperMACHO Project is an ongoing five-year microlensing survey of the LMC that 
is being carried out with the specific goal of answering the question, "Where do the lenses 
responsible for the measured event rates toward the LMC reside?" (Stubbs 1999). 

We have designed our survey to provide a significant increase in the number of detected 
events. This will allow a more robust assessment of the spatial variation of the microlensing 
optical depth across the face of the LMC, and will clarify whether the observed optical 
depth can be accounted for by LMC self-lensing — the most popular alternative to lensing 
by MACHO's. This paper presents an appraisal of SuperMACHO's ability to accomplish 
this goal. This assessment is based on extensive simulations that use observational data 
obtained during the first observing season. The LMC luminosity function plays a prominent 
role in this calculation, and we present an extensive analysis of this in a future paper (Rest 
et al. 2006). 

1.1. A First Step: Foreground lenses or LMC lenses? 

As a first step toward determining the nature of the lensing population, we will consider 
and evaluate two lensing scenarios: 
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• Screen-lensing : Lensing caused by a uniform (on the angular scale of the LMC) 
foreground lensing population. Examples are lensing by the Galactic halo, thick disk, 
or a many degree scale intervening population of lenses. 

• LMC Self-lensing : Lensing where the lens population is either the same as the 
source population or spatially close to the source population. Examples are LMC disk- 
disk, disk-bar, and bar-bar lensing and lensing of the LMC disk by a tidal tail within 
the LMC. 

We consider it a sensible first step to ascertain the extent to which these may be responsible 
for the microlensing events seen toward the LMC. 

The lensing rate for sclf-lensing shows a strong spatial dependence (e.g. the lensing 
rate for LMC bar-bar lensing is proportional to A''^^^, where is the areal density of stars), 
whereas the lensing rate of screen-lensing is directly proportional to the number of source 
stars observed. The goal of the SuperMACHO project is to determine which or what mixture 
of these two categories causes the observed microlensing rate. One key ingredient to achieving 
this is to increase the number of detected events. We do this in two ways: 1) by increasing 
both the number of source stars monitored, and 2) increasing our event detection efficiencies 
by performing difference image analysis. The corresponding improvement in event detection 
rate should move us out of the realm of small number statistics and allow us to determine 
the spatial distribution of the events on the sky (see Figure 1). This, in turn, should allow 
us to discriminate between the two possibilities described above. 

Our approach in assessing the survey's discrimination capability is to: 

1. Use actual LMC images obtained with the survey instrumentation to obtain star count 
information for each of our fields; 

2. Add simulated microlensed flux to the frames and assess the survey's event detection 
efficiency for each field as a function of input event parameters; 

3. Estimate, for the observed LMC optical depth, the likely event distribution statistics 
across the different fields under different lensing scenarios; 

4. Given the anticipated event detection rates, devise statistics that maximize the survey's 
ability to discriminate between screen and LMC self-lensing; and 



5. 



Assess the SuperMACHO survey's sensitivity to a specific illustrative LMC self-lensing 
scenario, namely the displaced LMC bar model proposed by Zhao & Evans (2000). 
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These steps are laid out below. In § 3 we describe the SuperMACHO survey strategy 
and the image analysis pipeline of the project. § 4 summarizes our parameterization of the 
LMC stellar luminosity function, an essential ingredient needed to go from the observable 
quantity, the number of microlensing events detected, to the physical properties of interest, 
in particular the mass and spatial distributions of the lens population. We present detailed 
luminosity function analysis in a companion paper. Rest et al. (2006). In § 5, we then use 
these results to model the number of detected microlensing events for the different candidate 
populations, and we predict whether it will be possible for SuperMACHO to distinguish 
between screen- and self-lensing. Before we begin a detailed examination of SuperMACHO's 
discrimination capability, we begin by examining the various lensing scenarios in greater 
detail. 



2. Lens Population Candidates 

The known and expected components of the Galactic/LMC system each contribute to 
the event rate toward the LMC. In the following we will discuss the different lens population 
candidates, each of which is problematic in some respect. 

2.1. Milky Way Halo Lenses? 

If the lenses reside as MACHO's in the Galactic halo, their inferred typical mass is 
0.1 — IMq (Alcock et al. 2000). If we assume that such a lens population is comprised 
of some known astronomical object, the most likely candidates are white dwarfs. There 
are some indications that there might be a previously undetected population of old white 
dwarfs (WDs) in the Galaxy (Ibata et al. 1999, 2000; Mcndez & Minniti 2000; Oppcnheimer 
et al. 2001; Nelson et al. 2002) favoring this interpretation. However, Kilic et al. (2004) has 
shown that some of these WD candidates are background AGN. Like any scenario populating 
the Galactic halo with stellar remnants, the halo WD explanation is challenged by stellar 
formation and evolution theory: the stellar progenitors are expected to enrich the gas and/or 
stars to a greater degree than has been observed. In addition, if other galaxies have similar 
halos, then their progenitor populations should be observable in galaxies at high redshift. 
There are ways out of these constraints, e.g. by assuming non-standard initial mass functions 
(Chabrier ct al. 1996; Chabrier 1999), or by anticipating lower metal yields from old, low- 
metallicity main sequence progenitors, or perhaps by allowing that the processed cjecta 
remain in the form of hot gas as yet undetected (Fields ct al. 2000). All these attempts require 
fine tuning of the models or invoke unlikely physics, rendering them somewhat unsatisfying. 
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The interpretation that the detected faint WDs are members of the Galactic halo is certainly 
not uncontested (e.g., see Richer 2001, Creze et al. 2004). 

2.2. Thick Disk Lenses? 

An alternate interpretation is that the lenses belong kinematically to the thick disk (Reid 
et al. 2001). In this scenario, the inferred number of WDs exceeds the number expected from 
the known thick disk or spheroid, forcing the invocation of an undetected very thick disk 

as an alternative to a halo population of lenses (Gates et al. 1998). Gyuk & Gates (1999) 
showed that such disks may be able to reproduce the observed optical depth toward the LMC. 
More recently, they showed that the predicted properties of such a population are consistent 
with the observed properties of the WDs (Gates & Gyuk 2001). Since the total mass of such 
thick disk WDs needed to account for the observed optical depth seen toward the LMC is 
much smaller than the total mass needed in the halo, this explanation solves some of the 
stellar evolution and chemical enrichment problems. Only more detailed observations can 
determine to which populations the WDs might belong or if they reside in a new, unknown, 
component of the Galaxy. 

2.3. LMC Self-lensing? 

A third possible interpretation of the optical depth toward the LMC is that the lenses 
are not part of our Galaxy but rather of the LMC itself, i.e. lens and source population re- 
side within the LMC. This was first suggested independently by Sahu (1994) and Wu (1994) 
and is denoted as LMC self-lensing. The most common self-lensing scenarios invoke pairs 
of LMC bar, disk, and halo objects as source and/or lens populations. Several groups find 
self-lensing optical depths close to 10^^ and claim that therefore the observed optical depth 
can be explained with self-lensing (e.g. Aubourg et al. 1999; Zhao & Evans 2000). That 
claim has been disputed by several other groups (e.g. Sahu 1994; Gould 1995; Alcock et al. 
1997; Gyuk et al. 2000) who find optical depths closer to 10~^ and thus not sufficient to 
explain the observed optical depth. The main differences between the different estimates of 
the optical depth come from different choices of LMC models and model parameters (Gyuk 
et al. 2000). Recent observations indicate that there might be kinematic (carbon star sample, 
Graff et al. 2000; RR Lyrae stars, Minniti et al. 2003) and photometric subcomponents of 
the LMC (Weinberg & Nikolaev 2001; van der Marel & Cioni 2001; van der Marel 2001). 
The unvirialized subcomponents can be caused by the tidal interaction between the LMC, 
the Galaxy, and/or the SMC. Some theoretical models that invoke such unvirialized sub- 
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components find that the optical depth is significantly increased and may account for half or 
even all of the microlensing event rate (Graff et al. 2000; Zhao & Evans 2000). The predicted 
events show peculiarities in their photometric, kinematic and spatial distribution, which can 
be used to distinguish between LMC stars and the other lens population candidates. For 
example, one of the two near-clump MACHO events (MACHO-LMC-1) is a few tenths of 
a magnitude fainter than the clump, and Zhao et al. (2000) argue that this is suggestive of 
having the lensed source star in a distinct population spatially separated and behind the 
LMC (and possibly more reddened). Using Hubble Space Telescope (HST) observations, Al- 
cock et al. (2001a) do not find any significant evidence for such systematically redder source 
stars. Their results marginally favor halo lensing. In addition, recent observations do not 
show any significant signs of kinematic outliers in the LMC (Zhao et al. 2003), restricting 
any additional kinematically distinct population to less than the 1% of the LMC stars. 



2.4. Galactic Halo Substructure? 

There has been increasing evidence that the Galactic stellar halo is not smooth. Beside 
the Magellanic Stream, another full-fledged tidal stream, the Sagittarius Dwarf galaxy, which 

currently is passing through the Galactic disk, has been detected (Ibata et al. 1995; Yanny 
et al. 2000; Ivezic et al. 2000; Vivas et al. 2001; Ibata et al. 2001). There is also tentative 
evidence for other tidal streams in the Galactic Halo (Newberg et al. 2002) and between the 
Galaxy and the LMC (Zaritsky & Lin 1997). Such an intervening dwarf galaxy or tidal tails 
could also cause the high microlensing event rate (Zaritsky & Lin 1997; Zhao 1998; Weinberg 
& Nikolaev 2001). 



3. Observations: The Implementation of the SuperMACHO Survey 

The primary motivation for starting the SuperMACHO survey was the collection of a 
sufficient number of microlensing events toward the LMC so that a statistical analysis can 
be made of competing theories for the location of the lenses. Previous LMC microlensing 
surveys such as MACHO and EROS have showed that microlensing can be detected and 
characterized with photometric sampling every few days and that LMC microlensing events 
are not shorter than two weeks. These surveys also highhghted the benefits of good seeing 
to reduce the effects of blending of source stars, while the simultaneous collection of data 
in multiple pass-bands to assess the achromaticity of candidate events has not proved very 
useful due to the effects of blending. With these lessons in mind, the SuperMACHO project 
was proposed to use a larger aperture to detect fainter events, at a better seeing site, in 



-9- 



a single filter and fit into the restrictions of using a non-dedicated telescope by observing 
every second night. The image analysis was designed to use difference image photometry 
which had been shown by previous surveys to be more efficient in detecting microlensing. 
The SuperMACHO proposal ^ was allocated 150-half-nights, distributed over 5 years, on the 
Cerro Tololo Interamerican Observatory (CTIO) Blanco 4m telescope through the NOAO 
Survey Program. The survey started in 2001 and will run through 2005. We note that we 
have waived any proprietary data access rights and that the SuperMACHO survey images 
are accessible through the NOAO Science Archive on the NOAO web site^. 

Observations arc carried out every other night in dark time during the months of Oc- 
tober, November, and December, when the LMC is most accessible from CTIO. We use the 
8K X 8K MOSAIC II CCD imager with a FOV of 0.33 square degree. The 8 SITe 2K x 4X 
CCDs are read out in dual-amphfier mode (i.e. different halves of each CCD are read out in 
parallel through separate amplifiers) to increase our observing efficiency. In order to max- 
imize the throughput we use a custom-made broadband filter (VR filter) from 500nm to 
750nm. The atmospheric dispersion corrector on the MOSAIC 11 imager allows for the use 
of this broad band without a commensurate Point-Spread-Function (PSF) degradation. 

In devising an observing strategy we want to find a good balance between maximizing 
the number of events detected and assuring a uniform spatial coverage. The work described 
here has guided our decisions on how to best spend the telescope time on the sky. We 
have defined a grid of 68 fields over the face of the LMC. Previous microlensing surveys 
found that the distribution of microlensing event durations^ toward the LMC has its peak 
at about i = 80 days, with virtually no event lasting less than two weeks. In order to sample 
the lightcurves adequately and sufficiently, we observe all fields every other night during 
dark time. This also serves to equafize, to first order, the event detection efficiency due to 
sampfing effects across the fields. There remains the field-dependent detection efficiency due 
to the different stellar densities and due to intentional inequality in exposure times. 

The distribution of the available observing time in a half-night across the LMC is driven 
by two confiicting considerations: the need to maximize the number of stars that we monitor 
and the need to survey as large a region as possible in order to discriminate between the 
different candidate lens populations. If the only goal were to maximize the number of 
monitored stars and, consequently, the number of detected microlensing events, we would 
concentrate on deep exposures of the central region of the LMC. Maximizing the number 



^http: / / www.ctio.noao.edu/supermacho 
^ftp: / / archive.tuc.noao.edu/pub/ 

^The event duration is twice the Einstein radius crossing time 
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of microlensing events is, however, less important than achieving maximum discrimination 
between models. 

We have adopted the strategy outlined by Gould (1999) to achieve a distribution of 
exposure times that maximizes the microlensing event rate subject to the constraints of a 
given spatial sampling and the differing sensitivity between the inner and outer regions of 
the LMC. The basic idea is that the distribution of exposure times for a microlensing survey 
is optimized when a shift of 5t in exposure from field A to field B gains as many stars in B 
as arc lost in field A. At this extremum = This condition must be achieved subject 
to two constraints. First, the total exposure time plus the total time spent on readout must 
equal the number of workable hours in a half night. The second constraint is that the number 
of sources monitored in the inner and outer regions of the LMC must be balanced to achieve 
the desired spatial coverage. 

We have used the stellar density normalizations described in § 4 and a simple division 
into inner and outer fields to optimize the distribution of exposure times given the properties 
of the MOSAIC II imager on the CTIO 4m telescope. We have set a minimum exposure 
time of 25 seconds in order to assure coverage of the sparser fields and a maximum exposure 
time of 200 seconds in order to avoid saturation effects. 

The Super MACHO Project started observations in Sept 2001. The data analysis pipehne 
is currently implemented as a combination of C code, IRAF, Perl and Python scripting tied 
together to provide an integrated but modular environment (Smith et al. 2002). 

The first steps of the data processing, crosstalk correction and astromctric calibration, 
are best done on the whole image because the units are not completely independent. The 
rest of the image reduction, as well as all of the transient analysis, breaks down naturally 
into 16 independent units, the amplifier-images^, and can therefore be efficiently handled in 
parallel. We employ a cluster of 18 CPUs with a 6.5 Terabyte redundant disk array. 

Standard photometry of transient or variable objects becomes inefficient in highly crowded 
images; therefore, we use a method called difference image analysis (DIA) which has rapidly 
evolved in the last few years. The first implementation was by PhiUips & Davis (1995) who 

introduced a method that registered images, matched the point spread function (PSF), and 
matched the fiux of objects in order to detect transients. Derivatives of DIA have been 
widely applied in various projects (e.g., MACHO, Alcock et al. 1999a; M31 microlensing, 
Crotts et al. 1999; OGLE, Wozniak 2000; WcCAPP, Gossl & Riffeser 2002; DLS, Becker 
et al. 2004). Since the PSF varies over the field-of-view due to optical distortions or out-of- 



''The MOSAIC II has 8 CCDs. Each CCD is read out by two amphfiers 
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focus images, for example, it is essential to use a spatially varying kernel (Alard & Lupton 
1998; Alard 2000). 

One of the main problems with the image differencing approach is that there are more 
residuals, e.g. cosmic rays and bleeds, than genuinely variable objects in the difference 
image. Therefore a standard profile-fitting software like DoPHOT (Schechter et al. 1993) 
has problems determining the proper PSF used to perform photometry in the difference 
image. When the difference image is analyzed with our customized version of DoPHOT, we 
force the PSF to be the one determined for the original, flattened image. Applying this a 
priori knowledge of the PSF helps to guard against bright false positives, such as cosmic rays 
and noise peaks, which generally do not have a stellar PSF. 

All detections are added into a database. Once the database is loaded and objects have 
been identified, queries are performed on new objects which are then classified. Objects of 
interest are then passed to a graphical user interface displaying stamps from image, template, 
and difference image for visual classification and interpretation. 

4. LMC luminosity functions 

In order to determine the optical depth from the observable quantities (the number 
and duration of detected events), the number of potential source stars must be known. 
This depends on both exposure depth and on the luminosity function (LF) of the source 
star population. Knowing the LMC LF is essential for the analysis (and prediction) of 
microlensing event rates. 

For this analysis we require the true LF for the stellar population of the whole LMC. The 
MACHO survey used Hubble Space Telescope (HST) images (Alcock et al. 1999b, 2001b) to 
determine the LF down to F ~ 24 for selected bar fields. They found that for all fields the 
LF was well-fit with a power-law with identical slope for V < 22.5, leveling out for fainter 
magnitudes. Even the HST LF shows a spread at magnitudes fainter than 22.5 (see Figure 
2 and 3 of Alcock et al. (2001b)). Since our survey is most sensitive to microlensing events 
with source star magnitudes in the range of 22 to 25 (see upper panel of Figure 6), these 
differences are important. Instead of relying upon a single LF for our analysis, we decided 
to explore a variety of LP's so that we can quantify the impact of choosing an incorrect LF 
on our conclusions. We used five different LP's to represent the possible range of LP's in 
the LMC or to represent possible variations in the global LF. Two of these LP's are the 
"limiting cases", while the three intermediate LP's are based on either LP's from the solar 
neighborhood or direct fits to our LMC photometry. As explained below, all five LP's are 



tied to a single power- law fit to the bright end of the LF, and the faint-end span a plausible 
range of luminosity distributions. We present our detailed analysis of the LMC luminosity 
function in a companion paper (Rest et al. 2006) and summarize its results here. 

We divided each of the 68 fields that we observe into 16 subfields based on the area 
covered by the MOSAIC II amplifiers (see Figure 1). For each of these subfields, we deter- 
mined an independent LF. First, we fit a single power-law with slope f3 and stellar density 
parameter $o (stars per square arcmin) to the bright end of the LF with a superimposed 
Gaussian function representing the red clump: 



where Njic, M^c^ and oj^c ^'^^ the total number, the average magnitude, and the spread of 
red clump stars, respectively. The instrumental magnitude M is defined as 



where / is the fiux in counts. Figure 2 shows the observed LF (green circles) for the mag- 
nitude range in which the completeness is close to 100% for amplifier 4 of field sme9. The 
dotted line (denoted as $>) is a completeness-corrected fit of $*(M) to the data in the 
magnitude range Mjic — 2.5 < M < Mjic -1-2.5. For fainter magnitudes (M > —9.5) 

the fit clearly overestimates the number of observed stars. This break in power-law is well 
known and documented for the solar neighborhood (Kroupa ct al. 1993; Rcid & Hawley 
2000; Chabrier 2001), and it is also seen in the LMC (Holtzman et al. 1997; Hunter et al. 
1997; Olsen 1999; Alcock et al. 2001b). Thus fitting a single power-law $> to the bright end 
of the LF's sets the stellar density parameter for each subfield and is the upper bound to 
plausible LF's. 

Since the faint end deviates noticeably from the single power-law, we parametrize the 
LF's as a combination of three power-laws: 



where M is defined in Equation 2. Some of our candidate LF's are only single or double 
power-laws. In these cases we disregard the functions which do not apply. We use the fitted 
red clump peak magnitude as the anchor point for the break magnitudes, i.e. Mi and M2 
are always with respect to Mrc- 

The following is a description of how we construct the five LF's: 




(1) 



M = -2.5 log/ 



(2) 




M < Ml 

Ml < M < M2 

M > M2 



(3) 
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• Upper limit LF ($>): The upper limit of the LF (black, dotted hne in Figure 2) is 

set by assuming that the power-law $*(Af) (sec Equation 1) fitted to the bright end 
continues to the faint end without a break ($i = ^q, Pi — (3 for all M). 

• Local neighborhood LF ($a) semi- empirical mass-luminosity relation: Reid 
et al. (2002) estimates a three power-law function as one of the possibihties of the 
present-day mass function using a semi-empirical mass-luminosity relation. We con- 
vert this present-day mass function back into a LF using the mass-luminosity func- 
tion for lower (Delfosse et al. 2000; Reid et al. 2002) and brighter (Reid et al. 2002) 
main-sequence stars, giving f3i = 0.34, j32 = 0.16, and (3^ = 0.27 with the two break 
magnitudes at V = 3.91 and V = 7.11. Using an unreddened absolute V magnitude 
of the red clump peak of 0.39 (Olsen, private communication), we can express the 
break magnitudes relative to the red clump as Mi = M^c + 4.3 and M2 — Mrc + 7.5. 
For each subfield we determine the stellar density parameter ^ao and the red clump 
peak magnitude Mrc by fitting the single power- law $*(M) (see Equation 1) to the 
bright end of the LF with the fixed slope oi P — Pi — 0.34. This LF is shown as red, 
short-dashed line line in Figure 2. 

• Local neighborhood LF empirical mass-luminosity relation: Another 
proposed form of the present-day mass function by Reid et al. (2002) is a two power- 
law based on an empirical mass-luminosity function. In the same way as for ^a, we 
estimate Pi — 0.38, P2 = 0.0314, and Mi = M^c + 5.0. The stellar density parameter 
$50 is determined by fitting the single power-law $*(M) to the bright end of the LF 
with the fixed slope oi P = Pi = 0.38. This LF is shown as red, long-dashed line in 
Figure 2. 

• Empirical LF {^c)'- As our empirical LF we fit the single power-law $*(M) to 
the bright end of the LF, yielding $1 = $0 and Pi = p. We fit a second power-law to 
the faint end of one of our sparse subfields, where the break in power-law is virtually 
unaffected by completeness (see Figure 2), yielding Mi = Mrc + 3.0 and P2 = 0.16. 
Note that the break 3 magnitudes fainter than the red clump peak magnitude is in 
very good agreement with the break seen in HST images of the LMC (see Alcock et al. 
2001b). Since we don't have data going deep enough to see the second break, we assume 
the same break magnitude and slope for the third power- law M2 — Mrc + 7.5 and 
Pz = 0.027, respectively, as determined for This LF is shown as blue, dot-dashed 
line in Figure 2. 

• Lower limit LF ($<): Similar to $> and $<7, we determine $1 and Pi by fitting the 
single power-law $*(M) to the bright end of the LF. As the first break in power-law, 
we use the same break we found for ^c'- Mi — Mrc + 3.0. As the slope we choose the 
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smallest slope of any of the power-laws, (32 — 0.027. This choice underestimates the 
LF at the faint end. This LF is shown as black, dot-dashed line in Figure 2. 

An overview of the different LP's is given in Table 1. In a future paper (Rest et al. 2006) we 
present a more detailed description on how we derive the LF's. 

We will use these different trial luminosity functions in the following section where 
we predict event rates in different microlensing scenarios. We show that a differential rate 
analysis of the data can discriminate between models independent of the actual underlying 
luminosity function. 



The SuperMACHO Project's initial goal is to distinguish between two broad categories 
of lensing: screen-lensing and self-lensing. We use data from the first year to predict the 
number of detected microlensing events from the different candidate populations in order 
to test if the SuperMACHO Project is able to distinguish between them. We add artificial 
stars to the images and determine the completeness of the detections. With that we estimate 
the spatial and temporal efficiency of detecting microlensing events in our data sets. In 
combination with the stellar luminosity functions described above, we predict and compare 
the number of microlensing events for various screen-lensing and self-lensing scenarios. 

The number of observed events depends on both the source and lens populations. In 
general, estimating the number of detectable microlensing events N^^i is complicated as it 
requires detailed knowledge about number density, velocity distribution, and other properties 
of the population. A good approximation is given by: 



where Nobs is the number of monitored stars for which microlensing can be detected, r is the 
microlensing optical depth, and £t is the sampling efficiency. This approximation separates 
the photometric and temporal completeness which makes the calculation much simpler. The 
number of monitored stars contains the photometric completeness and is 



where S{M) is the efficiency of detecting microlensing for a star with magnitude M, and 
$(M) is the luminosity function. The temporal completeness is contained in Ef as 



5. 



Event Rate Prediction 



Nmi = NobsrSt 



(4) 




(5) 




(6) 
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where T is the effective survey duration, i is the duration of a microlensing event, Prii) is 
the probabihty that the event is detected within T given the temporal cadence of the survey, 
and D{i) is a normahzed distribution of i. 

We calculate Nobs and £t independently for each field and amphfier, as described below. 

5.1. Number of Observed Stairs 

The best way to estimate the number of observed stars is to perform Monte Carlo 
simulations with all the images used. The shortcoming of this method is that it is very 
CPU-intensive and terabytes of images have to be simultaneously available on disk. Our 
goal is to predict a lower limit on the number of events; therefore, we choose a slightly 
different approach and perform Monte Carlo simulations on only a subset of images subject 
to certain assumptions. 

The question we have chosen to answer is: What is the probabihty £{M) that the flux 
of a star with a given magnitude M is magnifled during a microlensing event by a detectable 
amount? Let us assume that a microlensing event is detectable if the difference flux at peak 
has a signal-to- noise ratio S/N > 5. Whether such an event is then indeed detected depends 
on the temporal cadence as well as the seeing and transparency of the observing nights. We 
fold this into the temporal completeness analysis in § 5.2. 

During a microlensing event, the source star flux /o gets amplified to f = fo * A with 
an amplification that can be expressed as A{u) = [u^ + 2)/ {u\/u^ + 4), where u is the 
angular separation between source and lens in units of the Einstein angle. With difference 
imaging, the quantity we measure is the not the total amplified fiux, but the difference in 
ffux. A/ = f — fo, between the amplified and unamplified source star. By measuring only the 
difference flux, we avoid confusion due to crowding and blended sources. Using the difference 
images, then, we construct lightcurves of the difference flux. A/. 

Using Equation 2, we can then express this flux difference as an instrumental magnitude 
Mdiff where, 

Mdiff{M,u) = -2.51og(A/) (7) 
= M-2.51og(^(M) - 1) (8) 

We note that this is the instrumental magnitude of the flux difference, which is not the 
difference in magnitude between the amplifled and unamplifled source star. 

In § 6 we show that nearly all of the expected microlensing events have source stars 
in the magnitude range of 22.3 < Mvr < 25.3. Thus we can assume that the unamplifled 
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source star flux does not contribute significantly to the noise in the photometry, i.e. the 
photometry is background noise dominated. Within this limit, then, instead of adding the 
unamplified flux /q to the template image, and the amplified flux / to the image, we can 
just add the flux difference A/ to the detection image to test if we can detect such a flux 
difference. For this simulation we add flux into the images only and process these images 
through our difference image pipeline in the standard way. Because there is no source flux 
in the template, we measure the added difference flux. A/, in the resultant difference image. 
We then derive what fraction of the artiflcal stars we recover with a, S/N > 5 to obtain 
the empirical completeness function C{Mdiff). This allows us to estimate the probability 
that a microlensed star with intrinsic magnitude M has a change in flux whose S/N >5 ai 
maximum ampliflcation: 



where Uq is the angular separation at maximum brightness. With Equation 5, the number of 
observed stars Nobs can then be calculated using the luminosity function $(M) for a given 
held and amplifler (see § 4). 



The temporal cadence and observing conditions such as seeing and transparency have a 
signiflcant impact on event detection efficiency This is folded into P{t) and thus implicitly 
into St (see Equation 6). 

In order to estimate P{i) we perform Monte-Carlo simulations. During each year mi- 
crolensing events are drawn at random to have peak ampliflcation sometime during an inter- 
val of T = 300 days. This interval is dictated by the extent of the period of actual observing 
(100 days) and the need to ensure that real events which reach peak ampliflcation before 
or after that interval (but are observable during the survey months) are represented in the 
simulation. The 100-day padding at both ends of the observation period is driven by the 
characteristic i found by the MACHO Project. The SuperMACHO observation periods are 
usually allocated in three runs of about ten nights of bi-nightly observing separated by two 
weeks of bright time. Factors such as weather, instrument failures, and computer down time 
are taken into account by randomly eliminating one out of four nights. For a given i, we 
realize 1000 lightcurves for each microlensing impact parameter < uq < 0.5 in increments 
of 0.0125. By choosing uq — 0.5 as our upper limit we ensure that we have at least a 
magniflcation by a factor of 2.18. This will allow us to eliminate astrophysical sources of 
low-amplitude variability. Because the optical depth is typically calculated for Uq < 1, this 




(9) 



5.2. Sampling efficiency 
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simulation will only recover 50% of the number of events usually associated with a given 
optical depth The time of maximum is randomly chosen within T. As a lower hmit, we 
assume that all microlensing events observed have a difference flux with S/N {WQ,tQ) = 5, 
where Wq is the FWHM of the seeing at time to of maximum amplification. Then the S/N 
of a detection on another night at time t can be estimated as 



We randomly draw Wq and W from the distribution of seeing in the first year of the survey 
(see Figure 3). We want to emphasize that this is a conservative lower limit, since the stars 
we count in N^hg have a. S/N > 5 in the difference flux at maximum amplification; whereas, 
our Monte-Carlo simulations assume S/N — 5 at peak. 

For each hghtcurve we estimate the S/N ior every night we have taken data based on the 
FWHM and amphfication. The decision whether we do or do not detect such a hghtcurve 
is based on the following additional criteria: 

• At least 2 detections on the rising arm of the hghtcurve with a S/N > 2.0 

• At least 4 detections have a, S/N > 2.0 

This ensures that the events are "contained" within the survey coverage time and that there 
are enough significant detections to "trigger" an alert for follow-up observations. These 
cuts define the upper limit on the number of microlensing events we recover. We note the 
above set of criteria is insufficient to discriminate between microlensing and the population 
of background events in the actual survey. They are, however, comparable to the trigger 
criteria used in the MACHO alert system. Here we only consider the event sample in the case 
where there is no need for further discrimination. Discriminating microlensing is a separate 
problem and beyond the scope of this paper. 

In the final analysis, the initial cuts described above will define a population of candi- 
date events. We will apply additional cuts to eliminate contaminants such as supernovae, 
AGN's, and intrinsically variable stars. We will use photometric and spectroscopic follow-up 
observations to help define this advanced set of cuts. We note that these cuts will lower our 
detection efficiencies and require a reanalysis when they are determined. 

The upper panel of Figure 4 shows the probability P (i) of detecting a microlensing event 
with event duration of t. This event detection probability increases with event duration. Note 
that the probabihty is well below 1.0. This is because the interval used to define the input 
microlensing population was T — 300 days, much longer than the actual annual observing 



S/N {W, t) 




(10) 
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duration (~ 100 days), and the probability of detecting an event with a peak time to well 
outside the time the observations are taken is rather small. This effectively cancels out later 
on since we multiply by T when calculating St (see Equation 6). 

A more intuitive measure is P (i) x T/t, which is the number of microlensing events 
detected per l/r stars during T assuming that all microlensing events have an event duration 
of t (see lower panel of Figure 4). The decrease of detection probability for decreasing event 
duration is countered by that fact that for smaller event durations more microlensing events 
happen per l/r stars in the given observing time. This causes the peak at t ~ 50 days. 

Existing microlensing event statistics suggest that the event duration has a peak at 
about 80 days. Therefore we choose as D(t) a Gaussian distribution with the peak at 80 
days and with a spread of 20 days as inferred from Figure 9 in Alcock et al. (2000). Using 
Equation 6, we can now evaluate 

£t ~ 0.8 (11) 

Thus we will observe about 0.8 microlensing events per l/r observed stars in one survey year 
(spanning 100 days). 



5.3. Bar-Disk Self-lensing Models for the LMC 

The optical depth from self-lensing remains a matter of controversy. Some studies find 
rather small values in the range of (1.0 — 8.0) x 10~^ (e.g. Alcock et al. 1997, Gyuk et al. 
2000, and Jetzer et al. 2002); whereas, other studies suggest optical depths up to 1.5 x 10"^ 
(e.g. Zhao & Evans 2000). This controversy arises from the still rather imprecise knowledge 
of the structure of the LMC and consequent differences in the adopted models. 

The Zhao & Evans (2000) models derived in their paper are concrete and testable self- 
lensing models which we will use in the following sections to predict the SuperMACHO 
self-lensing event rate. As pointed out above, these models, denoted as model set A, see 
§ 5.3.1) predict a rather large optical depth and are thus more favorable to a self-lensing 
interpretation of LMC microlensing. We then improve upon their models (see § 5.3.2) for 
an alternative, more realistic model set, denoted as model set B. This allows us to make a 
direct spatial comparison between self-lensing and screen-lensing event rates. 



5.3.1. Zhao & Evans (2000) Model Set A 
Zhao & Evans (2000) derive in their paper concrete and testable self-lensing models. 
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They define a coordinate system with X, Y, and Z being decreasing right ascension, increas- 
ing decUnation, and hne-of-sight direction centered at the optical center of the LMC bar, 
where X, Y, and Z are in units^ of kpc. They approximate the average separation between 
source and lens as 

..^v-^ Ib^b + Ij^d + hh max( Afc + A^, A^^) 

where Id and are the star count density of the disk and bar, respectively. The line-of-sight 
depth of the bar and disk are denoted as A^ and A^, respectively, and A^^ is the hne-of-sight 
separation between the mid-planes of the bar and disk. In the limit in which the source and 
the lens are at roughly the same distance, the optical depth can then be expressed as 

^^^'^^ ^ '° 160 Me PC- Ikpc ^'^^ 

The value of r depends mainly on two parameters: the displacement Zq between the disk 
and the bar and the mass of the bar defined by the mass fraction /;,, which are implicit in 
the average separation and surface brightness. 

In § 6, we will vary these parameters to estimate the spatially varying optical depth for 
different model realizations, and we denote this set of models as model set A. 



5.3.2. Modified Zhao & Evans (2000) Model Set B 

The optical depth range for self-lensing found by Zhao & Evans (2000) is significantly 
larger than the ones found by other studies (e.g. Alcock et al. 1997, Gyuk et al. 2000, 
and Jetzer et al. 2002). As an alternative to model set A, we modify the original Zhao 
& Evans (2000) models by improving the effective separation A'(X, Y) (Nikolaev, private 
communication) which significantly decreases the optical depth: 

i^'^X.Y) = ^'A,/6 + /jA,/6 + hh^{^, + A,, A,,) ^^^^ 



where the function ^ is 



^i,(^ ^ \ \ \ = } ^bd ■ Afed > (Afe-F Arf)/2 

^ ' ^ A,, + ((A, + A,)/2-A,,)V(3A,A,) : A,, < (A, + A,)/^'^^ 



^For the LMC, 1 kpc is roughly 1 degree 
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We denote this model set as model set B. Note that the first two terms of A'{X,Y), the 
contributions of disk-disk and bar-bar self-lensing, arc smaller by a factor of six compared 
to the equivalent expression in the Zhao & Evans (2000) calculations (see Equation 12). 
This takes geometrical considerations into account since source and lens population are the 
same. The third term is a better approximation for the disk-bar separation: The two cases 
in Equation 15 represent situations 1) where disk and bar are well separated along the line of 
sight and 2) where they overlap each other. The two hmiting cases produce the same result 
when Abrf^ (Ab + Ad)/2. 

It is also important to note that the Zhao models have unrealistically high projected 
central surface densities, reaching 64OM0pc^^. A model with parameters as in Gyuk et al. 
(2000) leads to central densities on the order of 300 M© pc~^. The discrepancy is due to the 
very heavy bar (/;, = 0.5) and the quartic bar model (which is more centrally concentrated 
than a Gaussian bar) used by Zhao & Evans (2000). Less centrahzed projected surface 
densities result in a decrease of the central self-lensing optical depth. As before with model 
set A, we vary the displacement Zq and the mass fraction in § 6 and denote this set of 
models as model set B. 



6. Results 

In this section, we combine the results from the previous sections in order to obtain a 
quantitative prediction of the anticipated number of microlensing events. The main observ- 
able difference between self-lensing and screen-lensing is their distinctive spatial dependence 
of event rate; therefore, we calculate separately for each field and amplifier the expected 
number of observed microlensing events using Equation 4. 

First, we estimate the number of observed stars Nobs for each field and amplifier using 
Equation 5 as described in § 5.1. Figure 5 illustrates this for the example field sme9. The 
upper panel shows the range of luminosity functions we considered. The middle panel shows 
the detection efficiency^, and the lower panel shows the number of useful observed stars per 
magnitude, which is the product of the two upper panels. Two competing effects, increasing 
star counts but decreasing efficiency for fainter magnitudes, cause a peak at M ~ —7.6. Stars 
in the magnitude range — 9 < M < — 6 will contribute the most to the observed microlensing 
event rate. Using a zeropoint magnitude of 31.3 for this image, this corresponds to a VR 
filter (unlensed) magnitude range of 22.3 < Mvr < 25.3. The upper panel of Figure 6 shows 



^The efficiency leveis out at about 50% since we integrate in Equation 9 from to only 0.5, and not to 
1, due to tlie fact tliat we do not consider events witli amplification smaller than 2 
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the sum of the number of observed stars per magnitude bin for all fields and amplifiers, after 
correcting for exposure times. 

We obtain the total number of observed stars by integrating over magnitude (see lower 
panel). Not surprisingly, the LF chosen does significantly impact the estimated number of 
observed stars, e.g., choosing (short-dashed red line) produces ~ 2 times more observed 
stars than (dot-short-dashed blue line). The temporal completeness is taken into account 
by using £t ~ 0.8 (see § 5.2). 

We can now apply Equation 4 to calculate N^i for any combination of field, amplifier, 
and optical depth r. We use Tscreen — 1-2 x 10~^ for the optical depth for screen-lensing, 
as determined by the MACHO Project (Alcock et al. 2000). For sclf-lcnsing, we determine 
the optical depth for the Zhao & Evans self-lensing model sets A and B (see § 5.3) using 
Equation 13. The field-dependent Tse;/ for each of these models is calculated using each 
possible combination of /(, = [0.3,0.4,0.5], and Zq = [—1,-0.5,0,0.5,1], where /(, is the 
mass fraction of the bar, and Zq is the level of displacement between the disk and the bar in 
kpc. Based on observations, this covers the hkely parameter space of and Zq. The fields 
are divided into sets 1-5, based on their respective star density determined with the Zhao 
& Evans (2000) LMC bar model (see Figure 1). Set 1 is the most crowded (yellow) and 
set 5 is the least crowded (green). For each set of fields, we add up the predicted number 
of microlcnsing events for the different models and LF's (see Table 2). The upper panel 
of Figure 7 shows for each set of fields the number of microlcnsing events for Zhao LMC 
self-lensing model set A (black circles), model set B (blue squares), and screen- Icnsing (red 
triangle). For clarity the open symbols are plotted with a slight offset for a given set. Note 
that the spread in the anticipated number of detected microlcnsing events for screen-lensing 
is solely due to the different LF's ($^, and ^c) used (see Figure 2). For self-lensing, 
an additional source of spread is caused by using different values of and Zq. For sets 1-3 
the intrinsic difference in event rate for self- and screen-lensing is of the same order as the 
systematic errors. The event rate for self-lensing at the center of the bar is particularly large 
if the displacement between the disk and the bar is large {\Zo\ = 1). For the two outer sets, 
the rate of self-lensing strongly decreases and is well below the event-rate for screen-lensing. 



7. Discussion 

Generally, the intrinsic difficulty with drawing conclusions from the event rate is that 
there is a large spread in the predicted rate for a given field due to systematic biases, e.g. the 
LMC or Galaxy model used and the shape of the LF especially at the faint end. A way out 
of this dilemma is to go from an absolute measurement to a relative measurement: instead 
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of considering absolute event rates we investigate the differential event rate for a given field, 
defined as the ratio of event rate of the field to the total event rate of all fields (see lower 
panel of Figure 7). Using this normalized quantity to characterize the lensing rate across 
the LMC suppresses the dependence on luminosity function. This is clearly indicated in 
comparing the upper to the lower panels in Figure 7. To zeroth order the systematic error 
arising from LF uncertainty cancels out, and the measurement is much more robust. Table 3 
shows the differential event rates for self- and screen-lensing. 

Despite the anticipated increase in the number of microlcnsing events for the SuperMA- 
CHO survey, wc arc nevertheless faced with using small number statistics to try to distinguish 
between models for the lensing population. The basic approach we use is to consider what 
underlying rate could be statistically consistent with a given observed number of events. 

We represent these results in confidence level plots, for which we use the statistics 
given in Gehrels (1986). Because the differences between self- and screen-lensing are most 
pronounced in the outer field sets, we will investigate field set 5 and then the combined field 
sets 4 and 5. 

Figure 8 is an illustration of the outcome of this process. Assuming a total of 30 events 
are detected in the survey (our estimate of the lower bound we're likely to see for screen 
lensing), the x-axis corresponds to possible observed differential event rates in field set 5. 
The y-axis indicates the lowest allowable actual underlying differential rate, given Poisson 
statistics. The contours show the 90% (dot-dashed black line), 99% (dashed black hne), and 
99.5% (solid black line) confidence limits on the minimum actual underlying rate, given some 
measurement on the x-axis. The gray area shows the region excluded at 99.5% confidence. 

How can we use this plot? As a Gedankenexperiment, let us assume that screen-lensing 
is indeed the underlying mechanism. In that case, the expectation value for the field set 
5 differential event rate is 0.15 for the least favorable LF (see Column 3 in Table 3). Any 
particular experiment will measure a rate different than exactly 0.15, but it will most likely 
realize a value somewhere between 0.1 and 0.2. As a guide, this screen-lensing expectation 
value is indicated with a red dotted line in Figure 8. The question is, can these measurements 
exclude self-lensing? The upper limit of the differential rate of all Zhao self-lensing models in 
set A is 0.014 (see Column 5 in Table 3 and the blue, long-dashed line in the confidence plot). 
We can exclude these models with 99% and 99.5% confidence if we observe a differential rate 
larger than 0.11 and 0.13, respectively. Similarly for model set B, the maximum differential 
rate is 0.011 (blue short-dashed line), thus it can be excluded at 99.5% confidence if the 
observed differential event rate is bigger than 0.12. If SuperMACHO measures a field set 5 
differential rate larger than 0.13, than we can exclude self-lensing as the dominant mechanism 
at the 99.5% confidence level. 
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Let us consider the other case: If self-lensing is the dominant lensing process, then we 
expect to find many events in the central fields, and at most a couple in the outer fields. 
Therefore in this case it makes sense to calculate the upper bound of the allowed underlying 
differential rate given the observed number of events in the outer fields and in total. 

For example, the most optimistic self-lensing model A predicts 0.28 events in field set 5 

in 5 years. Most likely we will find cither or 1 event in this field set (assiiming self-lensing 
only). If we find more than 18 events in total none of which are in field set 5, then we can 
exclude screen-lensing (assuming a differential rate of 0.17, sec red, dashed line in the upper 
panel of Figure 9) at the 90% confidence level. Note that our sclf-lcnsing models predict a 
wide range of event rates from 5 to 35 events total during the 5 year survey (see Table 2), 
thus we can only expect to find that many self-lensing events if the LMC is similar to the 
models with strong central self-lensing. We would need 42 events to exclude screen-lensing 
at the 99.5% level and none in field set 5. 

If of 18 or more total events we observe none in field sets 4 and 5 (where screen lensing 
has an expectation rate of 0.30, see red, dashed line in middle panel of Figure 9), we can 
exclude screen-lensing with a confidence level greater than 99.5%. Finding no event is not 
unhkely since for the combined field set 4 and 5, the event rates predicted from our self- 
lensing models are between 0.25 and 1.6 events during the five years (see Table 2). Finding 
even one event in field sets 4 and 5 will weaken these conclusions, which can be seen by 
comparing the exclusion regions of the middle panel to those of the lower panel. 

In reality, the microlensing is probably not caused by a single lens population, but 

rather by a mixture of several populations. We can expect that it will be more difficult 
to differentiate between the populations. Still, very recent work re- investigating MACHO 
and ER0S2 events finds that even though some of the events are due to self-lensing, the 
total event rate and spatial distribution cannot be explained by self-lensing alone (Jetzer 
ct al. 2002). Also, a large spectroscopic survey targeting kinematic outhers in the LMC did 
not find evidence for a significant additional, kinematically distinct, population in the LMC 
(Zhao et al. 2003). The lack of such a population constrains the optical depth of self-lensing 
to values too small to explain the observed event rate. Even if screen-lensing is the cause 
for only a fraction of the observed event rate, we will be able to detect enough events in 
the outer field sets to exclude self-lensing as the sole lensing mechanism toward the LMC. 
If self-lensing is excluded as the sole lensing mechanism, the SuperMACHO event rate will 
provide a lower limit on the number of MACHO's in the halo of the MW and thus provide 
a lower limit of their contribution to the MWs dark matter. 



8. Summciry and Conclusions 



The reported microlensing event rate toward the LMC exceeds that expected from 
known visible components of our Galaxy, and the source of this observed excess rate is still 
the subject of discussion. Determining the nature of the lens population will have a great 
impact on our understanding of Galactic and LMG structure and possibly on the nature 
of dark matter. Possible explanations for the observed lensing can be broadly categorized 
into screen-lensing and self-lensing scenarios. Using the first year data of the SuperMAGHO 
Project, we performed completeness analysis by adding artificial stars to the images and 
estimated the spatial and temporal efficiency of detecting microlensing events in our data 
sets. We predicted a lower limit of observable microlensing events for both categories using 
the efficiency in combination with the stellar luminosity functions. We find that the Super- 
MAGHO Project should be able to distinguish between the two categories using the spatial 
differences in optical depth. Utihzing the differential event rate instead of the event rate 
itself decreases the impact of systematic errors rendering the results and conclusions more 
robust. 
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Fig. 1. — SuperMACHO fields superimposed onto the LMC. The fields are divided into sets 
1-5, based on their respective star density determined with the Zhao & Evans (2000) LMC 
bar model. Set 1 is the most crowded (yellow) and set 5 is the least crowded set (green). 
(LMC image courtesy of G. Bothun.) 
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Fig. 2. — This figure shows the different luminosity functions $>, ^a, ^b, ^Cj and $< for 
amplifier 4 of field sme9 versus the instrumental magnitude M. The units on the y axis 
are stars per square arcmin and magnitude bin. The black dotted line is the bright end 
LF function parametrized as a single power law ($>) fitted to the completeness corrected 
observed LF (green filled circles) at the bright end. 
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Fig. 3. — Seeing (FWHM) histogram for the 2001/2002 run for amphfier 4. The average is 
1.02" with a standard deviation of 0.14. 
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Fig. 4. — Probability of detecting a microlensing event with event duration of i. The upper 
panel shows the probability P (t) of detecting a microlensing event with event duration of 
i. The lower panel shows the more intuitive measure P (t) x T/i, which is the number of 
microlensing events detected per l/r stars during an interval T, assuming that all microlens- 
ing events have an event duration of i. The decrease of detection probability for decreasing 
event duration is countered by that fact that for shorter event durations more microlensing 
events happen per l/r stars in the given observing time. This causes the peak at f ~ 50 
days. 
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Fig. 5. — Determination of the number of observed stars in field sme9. The upper panel 
shows the LF candidates $a, ^b, ^c, and constructed as described in the text. The 
X-axis is the instrumental magnitude. The middle panel shows the event detection efficiency 
as a function of source star magnitude. The lower panel shows the number of monitored 
stars vs. magnitude, which is the product of the two upper panels. 
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Fig. 6. — Differential and cumulative Nobs for different LF's. The x-axis is the apparent VR 
magnitude, synthesized from standard V and R magnitudes through the equation VR = 
Ao - 2.5 X logio(v4ilO-°-^^ + (1 - Ai)10-°-^-^), where Aq and Ai are typically -1.1 and 0.37 
for the MOSAIC II, respectively. The upper panel shows the number of monitored stars per 
magnitude bin for all observed fields parameterized by luminosity function. The cumulative 
number of observed stars is then obtained by integrating over magnitude (see lower panel). 
The trial LF's $>, ^b, ^c, and $<, are described in the text. 
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Fig. 7. — Estimated microlensing event counts and differential event rates for the different 
field sets, as shown in Figure 1. The red triangles indicate screen-lensing for the different 
LF's. The Zhao & Evans LMC self-lensing model sets A and B are indicated with black 
circles and blue squares, respectively. The disk-bar displacement of the LMC models varies 
between -1 and 1 kpc, and the bar mass fraction varies between 0.3 and 0.5. For clarity the 
open symbols are plotted with a slight offset for a given set. The upper panel shows the 
number of events over a 5 year survey. The lower panel shows the differential event rate, 
normalized to the total number of events detected over the course of the survey. Note that 
the dependence on the LF's is greatly reduced for the differential rates, as shown in the lower 
panel. In particular, by comparing the innermost to the outermost fields we expect to be 
able to distinguish between the screen-lensing and LMC self-lensing scenarios. 
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Fig. 8. — Model exclusion plots for any observed differential microlensing event rate. The 
X-axis shows potential observed differential event rates for SuperMACHO field set 5. The 
y-axis corresponds to the allowed underlying event rate that SuperMACHO 's observations 
will use to constrain models. The upper limits to rates allowed by the observations at the 
90%, 99%, and 99.5% confidence levels are indicated with dot-dashed, dashed, and solid 
black lines, respectively, based on Poisson statistics assuming a total of 30 detected events. 
For a given observed differential rate (as plotted on the x-axis), the vertical projection of 
this measurement in the Figure will intersect the various confidence contours. The horizontal 
projection of this intersection to the y-axis yields the maximum differential ratio allowed by 
the observations. In particular, if the prior expectation value for screen lensing of 0.15 is 
actually observed by SuperMACHO, we can exclude with 99.5% confidence any model that 
predicts a ratio of 0.023 or less. This includes the self-lensing models A (0.014) and B (0.011) 
as described in the paper. 
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Fig. 9. — SuperMACHO's model discrimination ability as a function of the total number of 
detected microlensing events. The three panels correspond to different possible experimental 
outcomes (0 or 1 events in the outermost fields). The x-axis shows the total number of 
detected events. The y-axis shows the true differential event rate in fields 4 + 5, and in 
field 5 only. The lower boundary of the areas allowed by the observations at the 90%, 99%, 
and 99.5% confidence levels are indicated with the dotted, dashed, and solid black lines, 
respectively. For a given observed total number of microlensing events (x-axis) and the 
conditions listed inside each Figure window, the confidence contours show the maximum 
underlying differential rate allowed by the observations. For example, in the case where no 
event is observed in field sets 4 and 5 (middle panel), if 18 total events are observed, any 
model predicting an underlying differential event rate (set 4+5) greater than 0.3 (for example 
screen-lensing, see red, dashed line) is excluded at 99.5% confidence. 
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Table 1. Luminosity Function Parameterization 



LF 


$1 




/92 




Ml 


M2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


$> 


$0 


/3 














0.34 


0.16 


0.027 


+ 4.3 


Mrc + 7.5 


$B 




0.38 


0.031 




Mrc + 5.0 






$0 




0.16 


0.027 


Mrc + 3.0 


Mrc + 7.5 




$0 




0.027 




Mrc + 3.0 





Note. — Overview of the different luminosity function (LF) parameterizations used for 
computing the number of monitored sources, in the terminology used in the text. Column 
(1) shows the name of the LF model. Column (2) indicates the stellar density parameter 
used. Columns (3)-(5) indicate the relevant power law slopes, if applicable. Columns (6) 
and (7) show the transition magnitude between power-laws with respect to the fitted red 
clump magnitude Mrc- 
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Table 2. Predicted Microlensing Event Rates for the Super MACHO Survey 



Lens Pop. 


Field Set 


Min 


Mean 


Max 


(1) 


(2) 


(3) 


(4) 


(5) 


screen-lensing 


5 


4.2 


5.5 


6.9 




4+5 


8.2 


10.8 


13.8 




all 


25.2 


34.9 


45.2 


LMC self-lensing 


5 


0.12 


0.18 


0.28 


model set A 


4+5 


0.47 


0.87 


1.6 




all 


9.1 


18.7 


35.3 


LMC self-lensing 


5 


0.03 


0.07 


0.13 


model set B 


4+5 


0.25 


0.53 


1.13 




all 


5.2 


10.7 


21.4 



Note. — Predicted microlensing event totals for a 5 year SuperMACHO survey. The 
lens population and field set are indicated in column (1) and column (2), respectively. The 
minimum, mean, and maximum event rates for the different models and LP's are shown in 
column (3), (4), and (5), respectively. The Zhao & Evans LMC self-lensing model set A and 
B differ as described in § 5.3. Note that independent of the overall rate normahzation, the 
ratios of rates are a clear and robust indicator of the nature of the lensing population. 
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Tablc 3. Predicted Event Rate Ratios 



Lens Pop. 


Field Set 


Min 


Mean 


Max 


(1) 


(2) 


(3) 


(4) 


(5) 


screen-lensing 


5 


0.15 


0.16 


0.17 




4+5 


0.30 


0.31 


0.32 


LMC self-lensing 


5 


0.008 


0.010 


0.014 


model set A 


4+5 


0.038 


0.047 


0.055 


LMC self-lensing 


5 


0.004 


0.007 


0.011 


model set B 


4+5 


0.032 


0.051 


0.070 



Note. — Predicted differential microlensing event rates for screen and self-lensing. The 
lens population and field set is indicated in column (1) and column (2), respectively. The 
minimum, mean, and maximum differential rates for the different models and LP's are shown 
in column (3), (4), and (5), respectively. In order to be conservative, we use the minimum 
values for screen-lensing and the maximum values for self-lensing, indicated in bold. The 
Zhao & Evans LMC self-lensing model set A and B differ as described in § 5.3. Note that 
independent of the overall rate normalization, the differential rates are a clear and robust 
indicator of the nature of the lensing population. 



